function s=rick(f,nt,dt)
s=zeros(nt,1);
for i=1:nt
    tmp=pi*f*(i*dt-1./f);
    s(i)=(1-2*tmp^2)*exp(-(tmp)^2);    
end
end